Charge transfer statistics of a molecular quantum dot 
with a vibrational degree of freedom 
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We analyze the lull counting statistics (FCS) ol a single-site quantum dot coupled to a local Hol- 
stein phonon lor arbitrary transmission and weak electron-phonon coupling. We identify explicitly 
the contributions due to quasielastic and inelastic transport processes in the cumulant generating 
lunction and discuss their influence on the transport properties ol the dot. We find that in the 
low-energy sector, i.e. lor bias voltage and phonon frequency much smaller than the dot-electrode 
contact transparency, the inelastic term causes a sign change in the shot noise correction at certain 
universal values ol the transmission. Furthermore, we show that when the correction to the current 
due to inelastic processes vanishes, all the odd order cumulants vanish as well. 
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During the past decade molecular electronics has 
evolved into an important branch of condensed mat- 
ter physics Nowadays it is possible to electri- 
cally contact molecules of almost any geometry and 
complexity ranging from carbon nanotubes to hydro- 
gen molecules 0, H, 0, H, H, 01- Moreover, the experi- 
mental observables are not restricted to the linear con- 
ductance properties any more, but encompass also the 
nonlinear current-voltage characteristic I{V) as well as 
the shot noise Especially in the case when vibra- 
tional degrees of freedom are involved, these transport 
quantities display a number of very interesting features 

&l9|,lio|,[il[il,[il,[ii,[ii- 

One of the most pronounced effects of the interaction 
between the electronic and vibronic degrees of freedom 
is the abrupt change of the system conductance once the 
applied voltage is increased beyond a threshold which 
is related to the excitation energy of molecular vibra- 
tions i, i, i, [H E [S 113, H III- This is the reason 
why measurements around this turning point have be- 
come an invaluable instrument for the experimental in- 
vestigation of such systems. Interestingly, the sign of this 
conductance step depends crucially on the junction trans- 
parency, being negative for almost perfect transmission 
and positive in the opposite case. Most theoretical ef- 
forts have been centered around these two limiting cases 
@, H [13, [il S m Hi- According to Ref. [S], 
the transition is nonuniversal and the precise condition 
involves all system parameters 2J|. In general it occurs 
at an intermediate value of the transmission. 

Thus far mainly the nonlinear current-voltage charac- 
teristic has been analyzed in this regime. However, for 
future applications it is of importance to also possess 
information about the noise properties of such systems 
[1, 113] ■ A very convenient tool to calculate a variety of 
transport properties is the full counting statistics (FCS) 



which gives the probability distribution P{Q) to transfer 
Q elementary charges during a fixed (very long) waiting 
time T. The average value {{Q)) is then directly related 
to the current and its variance ((Q^)) to the noise power. 
The relevance of higher order cumulants ((Q")) has been 
demonstrated in Ref. [2^. Moreover, the analytic struc- 
ture of the cumulant generating function Inx(A) can give 
invaluable insights into the nature of t he p rocesses con- 
tributing to the transport [26, S H llSl • That is 
the reason why we would like to address the FCS of the 
molecular quantum dot coupled to a Holstein phonon. 
We model the system by the following Hamiltonian, 



H ^ Hr+HR + H, + HT + H, 



ph 



H, 



el.ph ■ 



(1) 



The terms Hl j^ describe the left and right electrodes 
in the language of the respective electron field operators 
We model them as noninteracting fermionic 
continua held at the chemical potentials ^l,r- The ap- 
plied bias voltage is then given hy V ~ ^ > ^ 
(we use units where e = h = 1). The particle transport 
between the dot and the electrodes is mediated by a lo- 
cal (symmetric |4l|) tunneling coupling at a; = with an 
amplitude 7 

Hr^ld) [M^ = 0) -I- ^r{x = 0)] + h.c. , (2) 

where d is the annihilation operator of the electron on 
the dot. We assume it to be spinless and modeled by 
Hd = Ad^d, A being the bare energy of the dot state 
which can be tuned by an applied gate voltage. The local 
(single- mode) phonon is described by Hph = fia^a while 
the electron-phonon coupling term with the amplitude g 
is given by 

Hei,ph = g q d^d, (3) 

where g = a + is proportional to the phonon displace- 
ment operator. 
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In order to determine the FCS, we calculate the cumu- 
lant generating function (CGF) Inx(A) = ln(e''^'^) which 
gives access to the cumulants (or irreducible moments) 
((Q")) of P{Q). The formalism for CGF calculation 
has been developed in Refs. [sil, [s^] and is by now well 
adapted to quantum impurity problems 3^, 33, 34|. We 
chose to use perturbation theory in the electron-phonon 
coupling g. The corresponding CGF then reads 



lnx(A)=lnxo(A)+lnx'(A). 

The first term is the CGF of the clean system at g 
At zero temperature, it is known to be given by 



(4) 



Inxo(A) = r ^ In [1 + Toiio^e^^ - 1) 



(5) 



where To(a;) = r'^/[{uj - A)^ + F^] is the single-particle 
transmission coefficient of the resonant level model |30|, 
35[ and F — 27r/9o7^ is the dot-electrode contact trans- 
parency, which in the wide flat band model depends on 
the energy independent density of electronic states in the 
electrodes po- The correction is treated in the spirit of 
Ref. [il] which is based on the generalized Keldysh ap- 
proach proposed in Ref. [1^. It is given by 



X'(A) 



c e 



-ig ds q{s) (s)d(s) 



(6) 



where the expectation value is taken with respect to the 
noninteracting Hamiltonian Hl + Hji + Hd + + Hph 
and is time-ordered on the Keldysh contour C. The de- 
pendence on the counting field A is contained in — 
^d) [e^^ipL + "^^Rj + li.c. The counting field has different 
signs on the two Keldysh branches, \{t) = ±A for t G dp. 



An application of the standard linked cluster expansion 
generates two different contributions to the lowest-order 
{g^) term: In^' = Inxi + ^^X2- The first one is present 
only if A 7^ 0, and is thus a consequence of a detuning of 
the dot from the particle- hole symmetric point. 



In XI (A) 



ig^T 



J2 ikl) dsi ds2 A'='(S1-S2) 



X N''is,)N\.S2). 



(7) 



Here, A{t,t') = —i{Tcq{t)q{t'))o denotes the phonon 
Green's function (GF) in Keldysh space and N''{s) = 
{d){tk)d{tk))\ - 1/2 for tk e Ck, k ^ ± are general- 
ized dot population probabilities. Due to the explicit 
A-dependence on the Keldysh branch k they are different 
on the forward/backward path Ck- 

The other contribution is given by 



lnX2(A) - E (fcO / ^7:'^\u;)A^'iu.) , (8) 

k,l=± •' 



where the generalized charge polarization loops are 
defined in terms of the dot level GFs D{t) = 
^i{Tcd(t)d)(t'))x by 



7r'='(a;) = I / ^D'^'iio + ij')D'\io') . (9) 



Both contributions ([7]) and® can be calculated using 
the unperturbed dot GFs [30| : 



D- 



D- 



(w - A) -t- iT{nL +nR-l) 



Voiio) V iT[e''^^^inL - 1) + e'^^^UR - 1)] ~{lu - A) + iTiuL + nR-l) 



(10) 



where the denominator is given by 

Voiiu) ^{lo- Af + F2 Ll(1 - nR){e'^ - 1) 



-I- nfl(l - nL)(e 



1) + 1 



(11) 



In the wide flat band model the electrodes are non- 
interacting fermions with Fermi distribution functions 
'^iJ,L(w). The chemical potentials in the left/right elec- 
trode are ±V/2. From now on we concentrate on the 
zero temperature results. The presented results hold 
for all temperatures much smaller than the smallest en- 
ergy scale in the system. The phonon GFs are given 
by A— (++)(i) = -iexp{Tifl\t\) and A-+^+-\t) = 



— z exp(±zrit). Equation PU)) allows one to calculate the 
generalized population probabilities N'' and one finds 
from Eq. ^ 



l-Xi(A) = -|^ 



In 



Ji-V/2)_ 
E.arctan(XZ^ 



(12) 



where we defined a A-dependent function f{uj) = [1 + 
To(ci;)(e*^ — 1)]~^. The contribution lnx2 can be written 
as a sum of two terms which describe quasielastic and 



inelastic processes [19']. The quasielastic part reads 



lnx«e;(A) 
where 



2tt 



duj- 

-y/2-A ^ 



2 4_ r2 e^-^ 



SgH, (13) 



kj=± 



(14) 



2r 



2fc //W2-A 
1 H arctan 



TT 



V(/F/2 - A)2 + r2 
+ (;i//2- A)| 



is the real part of the retarded dot self-energy S^(a;) [igj . 
It is related to the formation of the phonon sidepeak in 
the spectral function of the dot and can be considered as 
a renormalization of the bare transmission. 

The other contribution is the inelastic one[4^. It only 
contributes ioT: V > VL and is given by 



In Xmei(A) 



-e{v-n) I dw 



2tt 



y/2-A 



(15) 



E S{V/2 ~kA-ku!-n) 



1 



k=± 



w2 + r2e'^ 

iw(w + fcrj)-r2e*^"| 



uj{lu + fcil) 

(w + fcf7)2 +r2 ^ 2 (w + fcf7)2 + r2e*^ 



Equations p^ . and (fTS]) are the main result of this 
paper. Compatible results were obtained independently 
in two related studies 37 , [11] . They represent the full 
perturbative result for the correction to the CGF due to 
the phonon at arbitrary transmission and are valid at 
zero temperature and for small g (43j . 

The term In xi is due to the renormalization of the level 
energy A by the presence of the phonon 19], and is the 
leading term for small Q. This is precisely the condition 
when the adiabatic (or Born-Oppenheimer) approxima- 
tion is valid. The interpretation of other contributions is 
more lucid in the low energy sector when Q,V ^ T, A. 
Then we obtain 



2.g2TAy /AN 1 

^^XiiX) - -^arctan(-j ^^^^^^^. 



lnx,e/(A) = - 



g'rv 



1 



TT A2 -f r2 A2 + r2e*^ ' 



In Xtnel (A) 



g^T{V - n) 



TT 

A2 



A2 + r2 2 A2 + r2e 



e{v-n) 

1 A2 - r2e'^ 

2„iA 



, (16) 
(17) 
(18) 



A2 + r2e 



2„i\ 



In many cases the perturbative FCS of correlated systems 
is a linear combination of e*"^ terms with some integer n 
M, [13, m, m, [13]. One possible explanation is that 
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FIG. 1; (Color online) Density plot of the inelastic noise 
contribution S{{Q'^))inei/T for the parameters Q — Wg and 
V — 10. Ig as a function of F and A measured in units of g. 
Thin lines correspond to the transmissions given by Eq. (|21|l . 
Dark and bright areas represent negative and positive noise 
correction values, respectively. 



the corresponding terms describe transport processes in 
which the initial dot state is restored after the tunneling 
of n electrons. However, e*^ enters the above equations 
nonlinearly. This means that there are processes in which 
n electrons are necessary in order to bring the system 
back to its initial state: The dot can be excited by a 
first transmitted electron, then a number of them can 
fiow without interaction. The final nth electron then 
deexcites the system and leaves it in the initial ground 
state. 

An expansion of Eqs. (jlBUlSp in A/F, which corre- 
sponds to the limit of large transmission To{uj) sa 1, leads 
to a series containing e"™'*' {n > 0), describing elec- 
tron hacks cattering off the dot and thus a reduction of 
transmission. In the opposite limit of weak transmission, 
ro(a;) « 0, the CGFs can be expanded for small F/A. 
One then only encounters forward- scattering terms con- 
taining e™"^ [n > 0) which increase the transmission. 



Next we would like to discuss the issue of the sign 
change of the conductance and noise corrections in vicin- 
ity of = il. In the low-energy sector it is only governed 
by the inelastic term. The correction to the current is 
given by 



i d 
TdX 

9\y 



In Xmei(A) 



(19) 



A=0 

- r?) 



27r 



p2(p2 



(F2-hA2)3 



Here we immediately realize that the sign change occurs 
at precisely A = ±F. As the low-energy transmission 
coefficient is given by ro(0) = F2/(F2 -|- A2), the turning 
point condition is indeed To = 1/2. Interestingly, not 
only the current, but all the odd order cumulants (due to 
the inelastic correction) vanish at this point. 
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The correction to the noise power due to the inelastic 
tunnehng is given by 



5{{Q%nel 



A=0 



27r 



e{v-n) 



r2(r4_6r2A2 



(20) 



A* 




+ A2)4 

This term changes its sign at the values A = (±1 ± \/2)T 
which correspond to the bare transmission coefhcients 
given by 



(21) 



In Fig.[Tl we have plotted the inelastic noise contribution 
resulting from the CGF (fT5|) for a bias voltage V ^^l. As 
soon as r, A > y one clearly sees the change of sign. On 
the other hand, at low T,A<V, the observed features 
are nonuniversal. Since the sign change occurs only for 
V > and the elastic part is featureless around these 
critical values of A and F, we believe that the vanishing of 
the noise correction should be observable in experiments. 

This behavior proliferates to cumulants of higher or- 
ders. In the nth order the turning point condition is 
given by the zeros of an nth order polynomial. The 
corresponding solutions for the transmission coefficient 
are distributed symmetrically around Tq ~ 1/2. This is 
demonstrated in Fig. [2] One also notes that the cumu- 
lants increase drastically towards Tq ^ 1. Hence, these 
effects will be observable in junctions of the type used 
m Refs. 0,0,13, where large transmissions Tq « 1 have 
been observed. 

To summarize, we have calculated the charge trans- 
fer statistics of a molecular quantum dot with a single 
fermionic state coupled to a local Holstein phonon by 
means of a perturbative expansion in the electron-phonon 
coupling. We found that the FCS is the sum of an adia- 
batic (mean-field like) term, an elastic part and an inelas- 
tic term. The latter appears as soon as the applied volt- 
age exceeds the phonon frequency. We found that this 
inelastic term leaves a characteristic imprint on the noise 
power as well as on all higher order cumulants. We ex- 
pect these features to become observable in experiments 
in the nearest future. 
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